function xy=f_calibration_back(xyz,pars)
% 
CENTER_X=pars(1);
CENTER_Y=pars(2);
DX=pars(3);
DY=pars(4);
FOCALLENGTH=pars(5);
EP1=pars(6);
EP2=pars(7);
EQ1=pars(8);
EQ2=pars(9);

xy(:,1) = -FOCALLENGTH*xyz(:,1)./xyz(:,3)/DX+CENTER_X;
xy(:,2) = -FOCALLENGTH*xyz(:,2)./xyz(:,3)/DY+CENTER_Y;

end
